home *** CD-ROM | disk | FTP | other *** search
/ Aminet 15 / Aminet 15 - Nov 1996.iso / Aminet / misc / math / Laplace.lha / Laplace / Data / Projects.lha / LinA_Blatt14 < prev    next >
Encoding:
Text File  |  1996-08-26  |  1.7 KB  |  71 lines
  1. #LAPLACE
  2. #USEFLOAT=FALSE
  3. #USEERROR=FALSE
  4. #SIMPLIFY=TRUE
  5. #ERRBND=15
  6. #EXPDSP=5
  7. #PREC=8
  8. #F1=f(x)=
  9. #F2=g(x)=
  10. #F3=h(x)=f(x)/g(x)
  11. #F4=matrix(
  12. #F5=vector(
  13. #F6=diff(f)
  14. #F7=diff(g)
  15. #F8=diff(h)
  16. #F9=taylor(f,0)
  17. #F10=taylor(f,1)
  18. #F11=%(run type s:startup-sequence)
  19. #F12=
  20. #F13=
  21. #F14=
  22. #F15=
  23. #F16=
  24. #F17=
  25. #F18=
  26. #F19=
  27. #F20=
  28. m1=[1,2,-1,0;2,3,-2,2;-2,-4,0,5;0,1,2,3]
  29. m1b=invert(m1)
  30. m2=[2,-1,0,0;-1,2,-1,0;0,-1,2,-1;0,0,-1,2]
  31. const(a,b,c,d)
  32. m3=[0,-a,-b;a,0,-c;b,c,0]
  33. m=swaprows(m3,1,2)
  34. m=multrow(m,1,-inv(a))
  35. m=addrows(m,3,1,b)
  36. m=multrow(m,2,inv(a))
  37. m=addrows(m, 3, 2, c)
  38.  
  39. solve(m3)
  40. m4=[1,0,1,0;0,a,0,b;a,0,0,1;0,1,-b,0]
  41. det(m1)
  42. det(m2)
  43. det(m3)
  44. det(m4)
  45. m5=[a,-b,-c,-d;b,a,d,-c;c,-d,a,b;d,c,-b,a]
  46. d5=det(m5)
  47. rank(m1)
  48. solve(m1)
  49. rank(m2)
  50. solve(m2)
  51. rank(m3)
  52. solve(m3)
  53. rank(m4)
  54. solve(m4)
  55.  
  56. z1=[1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,1,0,0,0,0,0;0,0,1,0,0,0,0;0,0,0,1,0,0,0;0,0,0,0,1,0,0;0,0,0,0,0,1,0;0,0,0,0,0,0,1]
  57. z2=[1,0,0,0,0,0,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,0,1,0,0,0,0;0,0,0,1,0,0,0;0,0,0,0,1,0,0;0,0,0,0,0,1,0;0,0,0,0,0,0,1]
  58. z3=[1,0,0,0,0,0,0;0,1,0,0,0,0,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,0,0,1,0,0,0;0,0,0,0,1,0,0;0,0,0,0,0,1,0;0,0,0,0,0,0,1]
  59. z4=[1,0,0,0,0,0,0;0,1,0,0,0,0,0;0,0,1,0,0,0,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,0,0,0,1,0,0;0,0,0,0,0,1,0;0,0,0,0,0,0,1]
  60. z5=[1,0,0,0,0,0,0;0,1,0,0,0,0,0;0,0,1,0,0,0,0;0,0,0,1,0,0,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,0,0,0,0,1,0;0,0,0,0,0,0,1]
  61. z6=[1,0,0,0,0,0,0;0,1,0,0,0,0,0;0,0,1,0,0,0,0;0,0,0,1,0,0,0;0,0,0,0,1,0,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7;0,0,0,0,0,0,1]
  62. z7=[1,0,0,0,0,0,0;0,1,0,0,0,0,0;0,0,1,0,0,0,0;0,0,0,1,0,0,0;0,0,0,0,1,0,0;0,0,0,0,0,1,0;1/7,1/7,1/7,1/7,1/7,1/7,1/7]
  63. z_1=z7*z6*z5*z4*z3*z2*z1
  64. z_2=inv(z_1)
  65. z_2*[1,1,1,1,1,1,1]
  66. z_2*[7,6,5,4,3,2,1]
  67. z_3=combine(z_2-umatrix(7), [0,0,0,0,0,0,1])
  68. z_4=solve(z_3)
  69. cvector(z_4,7)
  70.  
  71.